A simple beam model for the shear failure of interfaces 
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We propose a novel model for the shear failure of a glued interface between two solid blocks. We 
model the interface as an array of elastic beams which experience stretching and bending under 
shear load and break if the two deformation modes exceed randomly distributed breaking thresh- 
olds. The two breaking modes can be independent or combined in the form of a von Mises type 
breaking criterion. Assuming global load sharing following the beam breaking, we obtain analyti- 
cally the macroscopic constitutive behavior of the system and describe the microscopic process of 
the progressive failure of the interface. We work out an efficient simulation technique which allows 
for the study of large systems. The limiting case of very localized interaction of surface elements is 
explored by computer simulations. 
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I. INTRODUCTION 

Solid blocks are often joined together by welding or 
glueing of the interfaces which are expected to sustain 
various types of external loads. When an elastic inter- 
face is subjected to an increasing load applied uniformly 
in the perpendicular direction, in the early stage of the 
failure process cracks nucleate randomly along the inter- 
face. Due to the heterogeneous microscopic properties of 
the glue, these cracks can remain stable under increas- 
ing load, which results in a progressive damage of the 
interface. This gradual softening process is followed by 
the localization of damage which leads then to the global 
failure of the interface and separation of the two solid 
blocks. 

Interfacial failure plays a crucial role in fiber reinforced 
composites, which are constructed by embedding fibers 
in a matrix material Q. Composites are often used as 
structural components since they have very good mass 
specific properties, i.e. they provide high strength with 
a relatively low mass, preserving this property even un- 
der extreme conditions. The mechanical performance 
of composites is mainly determined by the characteris- 
tic quantities of the constituents (fiber and matrix), and 
by the fabrication process which controls the material's 
microstructure, the formation of damage prior to appli- 
cations, and the properties of the fiber-matrix interface. 
In many cases the reinforcement is a unidirectional ar- 
rangement of long fibers resulting in highly anisotropic 
mechanical properties, i.e. in the direction of the fiber 
axis the composite exhibits high strength and fracture 
toughness since the load is mainly carried by fibers, how- 
ever, in the perpendicular direction the load bearing ca- 
pacity is provided solely by the matrix material. Hence 
the dominant failure mechanism of unidirectional com- 
posites perpendicular to the fibers' direction is shear. 
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Failure here occurs mainly due to the debonding of the 
fiber-matrix interface. 

Since disordered properties of the glue play a crucial 
role in the failure of interfaces, most of the theoretical 
studies in this field rely on discrete models 0, Q which 
are able to capture heterogeneities and can account for 
the complicated interaction of nucleated cracks. The pro- 
gressive failure of glued interfaces under a uniform load 
perpendicular to the interface has recently been studied 
by means of fiber bundle models 0, IE IE 0, IE 0, E3, EH - 
Several aspects of the failure process have been revealed 
such as the macroscopic constitutive behavior, the dis- 
tribution of avalanches of simultaneously failing glue and 
the structure of failed glue regions [12|]- Considering a 
hierarchical scheme for the load redistribution following 
fiber failure, a cascading mechanism was proposed for the 
softening interface in Ref . [l3L IT^ . The roughness of the 
crack front propagating between two rigid plates due to 
an opening load was studied in the framework of the fuse 
model. The microcrack nucleation ahead the main crack 
and the structure of the damaged zone were analyzed in 
detail . The shear failure of an interface between two 
rigid blocks has very recently been investigated by dis- 
cretizing the interface in terms of springs. It was shown 
that shear failure of the interface occurs as a first order 
phase transition [16|. 

In the present paper we study the shear failure of the 
glued interface connecting two solid blocks in the frame- 
work of a novel type of model. In our model the interface 
is discretized in terms of elastic beams which can be elon- 
gated and bent when exposed to shear load. Breaking of 
a beam is caused by two breaking modes, i.e. stretch- 
ing and bending, characterized by randomly distributed 
threshold values. The two breaking modes can be either 
independent or combined in terms of a von Mises type 
breaking criterion Assuming long range interaction 
among the beams, we obtained the full analytic solution 
of the model for the macroscopic response of the inter- 
face, and for the microscopic process of failure. We show 
that the presence of two breaking modes lowers the crit- 
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ical stress and strain of the material without changing 
the statistics of bursts of simultaneously failing elements 
with respect to the case of a single breaking mode. The 
coupling of breaking modes results in further reduction of 
the strength of the interface. We demonstrate that vary- 
ing the relative importance of the two breaking modes 
the macroscopic response of the interface can be tuned 
over a broad range. The limiting case of very localized in- 
teraction of beams is also considered. We determine the 
constitutive behavior and the distribution of avalanches 
of breaking beams for the case when beams interact solely 
with their nearest and next-nearest neighbors in a square 
lattice. An effective simulation technique is worked out 
which makes it possible to study systems of large size. 



II. PROPERTIES OF THE MODEL 

In our model we represent the glued interface of two 
solid blocks as an ensemble of parallel beams connecting 
the two surfaces. First, we derive an analytical descrip- 
tion of a single beam of quadratic cross section clamped 
at b oth e nds and sheared by an external force /, see 
Fig. |l(a)| The shearing is exerted in such a way that 
the distance I between the two clamping planes is kept 
constant. Consequently, the beam experiences not only 
a torque m, but also a normal force t due to the elonga- 
tion A/, which is characterized by the longitudinal strain 
e = Al/l. 
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FIG. 1: (a) Shearing of a single beam between two rigid plates. 
Since the distance I between the plates is kept constant, the 
beam experiences stretching and bending deformation, with 
longitudinal t and shear / forces, (b) Shearing of an array of 
beams, with the corresponding forces. In the case shown, one 
beam is broken. 

We derive the form of the deflection curve of the beam, 
as well as the magnitude of the tension force. It is nec- 
essary to introduce some approximations, so the model 
can be incorporated into the simulation code in a sensi- 
ble way. Following the procedure outlined e. g. in [l8| . 
we solve the differential equation for the beam deflection 



under the influence of the lateral force / and a given 
stretching force t. We then solve self-consistently for t(f), 
with t being the result of the longitudinal elongation. 

The governing differential equation for the bending sit- 
uation depicted in Fig. 1(a) can be cast in the form 
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EV 



(1) 



with boundary conditions 



C(o) = o, 
C'(o) =o, 

C"(l/2) =0. 



(2) 



Here, E denotes the modulus of elasticity, and I is the 
moment of inertia for bending of the beam. For a beam 
of rectangular cross-section, we have I = d 4 /12, where d 
is the side length. Let us briefly motivate this ansatz by 
stating that the second derivative C"i z ) is proportional to 
the torque on the beam, so consequently it needs to van- 
ish at the beam half-length 1/2. Accordingly, the third 
derivative C"'( z ) is proportional to the shearing force ex- 
erted on the beam, hence, it constitutes a term of the 
balance equation, Eq. Q • The first derivative term with 
C'(z) denotes the projection of the tension force t. Due to 
the clamping, the deflection and its first derivative must 
vanish at the end = 0. The formula for the bending 
moment m is 

m = -ElC'(z) . (3) 
The solution ((z) for vanishing t can be obtained as 0] 



CM = ^(31-2,), 
from which we can calculate the elongation 
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C 2 dz. 



It follows from the above equation 

t = ES^- = ESe, 



(4) 



(5) 
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where S = d 2 is the beam cross-section area. The first 
order solution for t(f) reads as 
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From a computational point of view, a formulation of 
bending and stretching in terms of the longitudinal strain 
e is more suitable than using the lateral force /. For that, 
we only need to replace m(f) by m(e), which yields 



m(e) 
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with 



3/ J 
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(9) 



Using e as an independent variable enables us to make 
comparisons to the simple case of fiber bundle models 
HIE H UmiH where the elements can have solely 
stretching deformation. In the model we represent the 
interface as an ensemble of parallel beams connecting the 
surface of two rigid blocks (see Fig. 1(b) I. The beams are 
assumed to have identical geometrical extensions (length 
I and side length d) and linearly elastic behavior charac- 
terized by the Young modulus E. In order to capture the 
failure of the interface in the model, the beams are as- 
sumed to break when their deformation exceeds a certain 
threshold value. As it has been shown above, under shear 
loading of the interface beams suffer stretching and bend- 
ing deformation resulting in two modes of breaking. The 
two breaking modes can be considered to be independent 
or combined in the form of a von Mises type breaking cri- 
terion. The strength of beams is characterized by the two 
threshold values of stretching ei and bending e 2 a beam 
can withstand. The breaking thresholds are assumed to 
be randomly distributed variables of the joint probabil- 
ity distribution (PDF) p(ei,e 2 ). The randomness of the 
breaking thresholds is supposed to represent the disorder 
of the interface material. 

After breaking of a beam the excess load has to be 
redistributed over the remaining intact elements. Cou- 
pling to the rigid plates ensures that all the beams have 
the same deformation giving rise to global load sharing, 
i.e. the load is equally shared by all the elements, stress 
concentration in the vicinity of failed beams cannot oc- 
cur. If one of the interfaces has a certain compliance, the 
load redistribution following breaking of beams becomes 
localized. This case has recently been studied for the ex- 
ternal load imposed perpendicular to the interface [l6| . 

In the present study we are mainly interested in the 
macroscopic response of the interface under shear load- 
ing and the process of progressive failure of interface ele- 
ments. The global load sharing of beams enables us to ob- 
tain closed analytic results for the constitutive behavior 
of the system for both independent and coupled breaking 
modes. We examine by computer simulations the statis- 
tics of simultaneously failing elements. The limiting case 
of the very localized interaction of interface elements is 
explored by computer simulations. 



III. CONSTITUTIVE BEHAVIOR 

Assuming global load sharing for the redistribution of 
load after the failure of beams, the most important char- 
acteristic quantities of the interface can be obtained in 
closed analytic form. 

Breaking of the beam is caused by two breaking modes, 
i.e. stretching and bending which can be either inde- 
pendent or coupled by an empirical breaking criterion. 



Assuming that the two breaking modes are independent, 
a beam breaks if either the longitudinal stress t or the 
bending moment m exceeds the corresponding break- 
ing threshold. Since the longitudinal stress t and the 
bending moment m acting on a beam can easily be ex- 
pressed as functions of the longitudinal deformation e, 
the breaking conditions can be formulated in a transpar- 
ent way in terms of e. To describe the relative impor- 
tance of the breaking modes, we assign to each beam 
two breaking thresholds e\, e\, i = 1, . . . , N, where N de- 
notes the number of beams. The threshold values t\ and 
62 are randomly distributed according to a joint proba- 
bility density function p(e\, 62) between lower and upper 
bounds e™" , e™ aa: and e™ m , e™ ax , respectively. The den- 
sity function needs to obey the normalization condition 



de 2 / dei p(e u e 2 ) = 1. (10) 



A. OR breaking rule 

First, we provide a general formulation of the failure 
of a bundle of beams. We allow for two independent 
breaking modes of a beam that are functions / and g of 
the longitudinal deformation e. Later on this case will 
be called the OR breaking rule. A single beam breaks if 
either its stretching or bending deformation exceed the 
respective breaking threshold e\ or £2, i.e. failure occurs 
if 



m 

9(e) 
£2 



> 1 or 

> 1, 



(11) 
(12) 



where Eqs. (|11I12I) describe the stretching and bending 
breaking modes, respectively. The functions /(e) and 
g(e) are called failure functions, for which the only re- 
striction is that they be monotonic functions of e. For 
our specific case of elastic beams the failure functions 
can be determined from Eqs. KilSI as 



/( e ) = e 1 5(e) = a\ 



(13) 



where a is a constant and the value of the Young modulus 
E is set to 1. 

In the plane of breaking thresholds each point (£1,62) 
represents a beam. For each value of e those beams which 
survived the externally imposed deformation are situated 
in the area /(e) < ei < e™ ax and g(e) < e 2 < ef ax , as 
it is illustrated in Fig. [2 Hence, the fraction of intact 
beams N- Ultact /N at a given value of e can be obtained by 
integrating the density function over the shaded area in 
Fig. El 



intact 

N 



de 2 / dei p(ei,e 2 ). (14) 
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FIG. 2: Plane of breaking thresholds (ei,e2). The point of 
intersection of /(e) and g(e) determines the fraction of re- 
maining beams. 



Due to the global load sharing, deformation and stress 
of the beams are the same everywhere along the inter- 
face. Consequently, the macroscopic elastic behavior of 
the system can be obtained by multiplying the load of a 
single beam, er^ = e (E = 1 is taken), by the fraction of 
intact elements Eq. (|14|l 



de 2 



9(0 



dei p(ei,e 2 ) 



(15) 



/(e) 



Assuming that the breaking thresholds, characterizing 
the relative importance of the two breaking modes, are 
independently distributed, the joint PDF can be factor- 
ized as 



£>(ei,e 2 ) =Pi(ei) -p 2 (e 2 ). 



(16) 



Introducing the cumulative distribution functions 
(CDFs) as 

Pi(ei) = f 1 pi(ei)dei, andP 2 (e 2 ) = f * p 2 (e' 2 )de' 2 , 

(17) 

we can rewrite Eq. I|15[) as 



(18) 



e I de 2 p 2 (e 2 ) / de 1 p 1 (e 1 ) 
= e[l-P 2 ( ff (e))][l-Px(/(e))]. 



This is the general formula for the constitutive behavior 
of a beam bundle with two breaking modes applying the 
OR criterion. In the constitutive equation 1 — Pi (/(e)) 
and 1 — P 2 (g(e)) are the fraction of those beams whose 
threshold value for bending and stretching is larger than 
17(e) and /(e), respectively. It follows from the structure 
of Eq. I|18|l that the existence of two breaking modes leads 
to a reduction of the strength of the material, both the 



critical stress and strain take smaller values compared to 
the case of a single breaking mode applied in simple fiber 
bundle models [ISBHSHIISIllI- 

Considering the special case of two uniform distri- 
butions for the breaking thresholds in the intervals 
j e mm £ maxj an( j j £ min^ £ max^ respectively, we can derive 

the specific form of Eq. I|18f) by noting that 



1 



1 



(19) 



After calculating the cumulative distributions, the final 
result follows as 



[£ max_ /(e)][e 



max 
2 



9(e)] 



max ,111111] [ max _minl 



(20) 



More specifically, if the distributions have equal bound- 
aries [0, 1], and substituting the failure functions / and g 
from Eq. (|13fl , the constitutive equation takes the form 



a= e[l - e][l - aVe\. 



B. Von Mises type breaking rule 



(21) 



We now address the more complicated case that the 
two breaking modes are coupled by a von Mises type 
breaking criterion: a single beam breaks if its strain e 
fulfills the condition [ljj 



f(e)Y , 5 (e) 



ci 



> 1. 



C-2 



(22) 



This algebraic condition can be geometrically represented 
as it is illustrated in Fig. [3] In the plane of the fail- 
ure thresholds ei, e 2 , the beams that survive a load e are 
bounded by the maximum values e5" ax , e™ 8 * and the hy- 
perbola defined by Eq. \2'2\) . Calculating the intersection 
points a and b defined in Fig. |3 which are found to be 



= m 



,max 
e 2 



9{e) 



1/2 



and 



9(e) (er x ) 2 



(£ max )2 _ /2(e ) ' 

the fraction of surviving beams can be expressed as 

N;„t„ r t 



(23) 



N 



dei I de 2 p(ei,e 2 ) (24) 

62(^1, e) 



with the integration limit 
e 2 (ei,e) = 



el9(e) 



(25) 



The constitutive behavior in this case is therefore given 
by 



dei I rfe 2 p(ei,e 2 ) . (26) 

'62(61, e) 
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FIG. 3: Intact beams in the plane of the failure thresholds 
ei,£2 for a given strain e, if breaking is determined by the 
von Mises criterion. The values a and b are defined as the 
intersections between the curve of the breaking condition Eq. 
lUD and the maximum values e" lax , e ™ ax : respectively. The 
shaded region labeled A denotes the intact beams; the shaded 
region B represents the additionally failing beams that would 
be intact in the case of the 0_R-criterion. 



FIG. 4: Constitutive behavior of a bundle of beams with two 
breaking modes in a strain-controlled simulation of N = 4- 10° 
beams, under the OR (dashed line), von Mises type (solid 
line), and a pure stretching breaking criterion (dotted line). 
The random failure thresholds for the breaking modes of each 
beam are sampled from uniform distribution between [0, 1]. 
The points marked with '+' and 'x' denote the respective 
theoretical results, Eqs. I21II27I . and a = e(l — e) for the pure 
stretching case. The constants E and a are set to unity here. 



We would like to emphasize that assuming independence 
of the breaking thresholds the joint distribution factor- 
izes p(ei,€2) = pi(ei) ■ ^2(^2), but the integrals in Eq. 
(|26[) over the two variables cannot be performed inde- 
pendently. Still, the integral in Eq. (|26|) can be evaluated 
analytically for a broad class of disorder distributions. As 
an example, we again consider two homogeneous distri- 
butions Eq. (|19|l over the interval [0, 1] along with the 
failure functions Eq. (|T3|) . Setting the Young modulus 
and the parameter E = 1 = a, the integrals yield 




(27) 



Even for the simplest case of uniformly distributed 
breaking thresholds, the constitutive equation takes a 
rather complex form. It is important to note that the 
coupling of the two breaking modes gives rise to a higher 
amount of broken beams compared to the OR criterion. 
In Fig. |21 the beams which break due to the coupling of 
the two breaking modes fall in the area labeled by B. 



In order to determine the behavior of the system for 
complicated disorder distributions and explore the micro- 
scopic failure process of the sheared interface, it is nec- 
essary to work out a computer simulation technique. In 
the model we consider an ensemble of N beams arranged 
on a square lattice. Two breaking thresholds e\,e% are 
assigned to each beam i (i = 1, . . . , N) of the bundle 
from the joint probability distribution p(e\, €2). For the 
OR breaking rule, the failure of a beam is caused either 
by stretching or bending depending on which one of the 
conditions Eqs. (|11I12|) is fulfilled at a lower value of the 
external load. This way an effective breaking threshold 
el can be defined for the beams as 



min(r 1 (ei)><T 1 (4)), 



i = l,. 



(28) 



where / _1 and g~ x denote the inverse of /, g, respec- 
tively. A beam i breaks during the loading process of the 
interface when the load on it exceeds its effective breaking 
threshold e l c . For the case of the von Mises type breaking 
criterion Eq. 11' 21) . the effective breaking threshold e* of 
beam i can be obtained as the solution of the algebraic 
equation 



9(4) 



1. 



1. 



(29) 
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FIG. 5: The beams that break due to mode g fall in the 
shaded region. The labels a and b mark the abscissae of the 
intersection points of the curve g~ 1 {e2) = /~ 1 (ei) with the 
lines £2 = e™" 1 an d £2 = £2^, respectively. 



FIG. 6: Constitutive behavior of a bundle of N = 90000 
beams using the OR criterion. The parameter values Ai = 1.0 
(stretching), mi = = 2 were fixed, while A2 corresponding 
to the bending mode was shifted. Inset: Fraction of beams 
breaking by stretching and bending as a function of A2. 



Although for the specific case of the functions /, g given 
by Eqs. (|11I12|) the above equation can be converted to 
a 4th order polynomial and solved analytically, this solu- 
tion turns out to be impractical, especially since the nu- 
merical evaluation of the solution is too slow. We there- 
fore solve Eq. I|29[) numerically by means of a modified 
Newton root finding scheme, where we make use of the 
fact that the solution has the lower bound 0. 

In the case of global load sharing, the load and de- 
formation of beams is everywhere the same along the 
interface, which implies that beams break in the increas- 
ing order of their effective breaking thresholds. In the 
simulation, after determining e* for each beam, they are 
sorted in increasing order. Quasi-static loading of the 
beam bundle is performed by increasing the external load 
to break only a single element. Due to the subsequent 
load redistribution on the intact beams, the failure of a 
beam may trigger an avalanche of breaking beams. This 
process has to be iterated until the avalanche stops, or 
it leads to catastrophic failure at the critical stress and 
strain. Under strain controlled loading conditions, how- 
ever, the load of the beams is always determined by their 
deformation so that there is no load redistribution and 
avalanche activity. 

In Fig. 0] the analytic results of Sec. IHII on the consti- 
tutive behavior Eqs. l|21l I27|) are compared to the corre- 
sponding results of computer simulations. As a reference, 
we also plotted the constitutive behavior of a bundle of 
fibers where the fibers fail solely due to simple stretch- 



ing H H II 0, II H G3 113 It can be seen in the figure 
that the simulation results are in perfect agreement with 
the analytical predictions. It is important to note that 
the presence of two breaking modes substantially reduces 
the critical stress a c and strain e c (a and e value of the 
maximum of the constitutive curves) with respect to the 
case when failure of elements occurs solely under stretch- 
ing. Since one of the failure functions g(e) is non- linear, 
the shape of the constitutive curve cr(e) also changes, es- 
pecially in the post-peak regime. The coupling of the 
two breaking modes in the form of the von Mises crite- 
rion gives rise to further reduction of the strength of the 
interface. 



V. PROGRESSIVE FAILURE OF THE 
INTERFACE 

During the quasi-static loading process of an interface, 
avalanches of simultaneously failing beams occur. Inside 
an avalanche, however, the beams can break under differ- 
ent breaking modes when the OR criterion is considered, 
or the breaking can be dominated by one of the breaking 
modes in the coupled case of the von Mises type crite- 
rion. Hence, it is an important question how the frac- 
tion of beams breaking due to a specific breaking mode 
(stretching or bending) varies during the course of load- 
ing of the interface. 
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FIG. 7: Constitutive behavior for different values of the shape 
parameter mi of stretching. Strain controlled simulation of 
N = 90000 beams with failure due to the OT?-criterion, fixing 
the parameters Ai = A2 = f.O and m2 = 2. Inset: total 
fraction of beams broken under mode g during the course of 
loading. 



For the OR criterion, those beams break, for instance, 
under bending, i.e. under mode g defined by Eq. I)12[l. 
whose effective breaking threshold is determined by 



a 



L (e 2 ) in Eq. i)28|) so that the inequality holds 



< 



(30) 



In the plane of breaking thresholds {£1,62} the region 
of beams which fulfill the above condition is indicated by 
shading in Fig. |3J The fraction of beams B g (e) breaking 
under mode g up to the macroscopically imposed defor- 
mation e can be obtained by integrating the probability 
distribution p(ei, e 2 ) over the shaded area in Fig. EI Tak- 
ing into account the fact that the intersection points a, b 
defined in Fig. [5] may in general lie outside the rectangle 
(e™ in ,e™ ax ,e™ ax ,e™ in ) and adjusting the integral limits 
accordingly, we arrive at the following formula for the 
fraction of fibers breaking under mode g as a function of 



min(/(e),6) ffCrHei)) 



B g {e)= J cfei J de 2 p(e 1: e 2 ) 

max(e'j mri ,a) e 2 im 
/(e) 

+ J dei J de 2 p(ei,e 2 ) 

min(/(e),b) e f in 
«™ aX 9(e) 

+ J dei J de 2 p(e 1 ,e 2 ). 



(31) 



It should be noted that the second integral vanishes un- 
less b < e™ ax . The total fraction of beams breaking under 
mode g during the entire course of the loading can be ob- 
tained by substituting e = e max in the above formulas, 
where e max denotes the deformation at the breaking of 
the last beam. 

In order to study the effect of the disorder distribu- 
tion p(ei,e 2 ) of beams on the relative importance of the 
two breaking modes and on the progressive failure of the 
interface, we considered independently distributed break- 
ing thresholds t\ , e 2 both with a Weibull distribution 



Pb{6b) = U 



exp 



(32) 



where index b can take values 1 and 2. The exponents 
mi,m 2 determine the amount of disorder in the system 
for stretching and bending, respectively, i.e. the width 
of the distributions Eq. (|3~2"|) . while the values of Ai,A 2 
set the average strength of beams for the two breaking 
modes. Computer simulations were performed in the 
framework of global load sharing by setting equal val- 
ues for the shape parameters mi = m 2 and fixing the 
value of Ai = 1 of the stretching mode, while varying A2 
of the bending mode. 

The total fraction of beams breaking by stretching and 
bending using the OR breaking rule is presented in Fig. 
El Increasing A2 of the bending mode, the beams become 
more resistant against bending so that the stretching 
mode starts to dominate the breaking of beams, which is 
indicated by the increasing fraction of stretching failure 
in the figure. In the limiting case of A2 >> Ai the beams 
solely break under stretching. Decreasing A2 has the op- 
posite effect, more and more beams fail due to bending, 
while the fraction of beams breaking by the stretching 
mode tends to zero. It is interesting to note that varying 
the relative importance of the two failure modes gives also 
rise to a change of the macroscopic constitutive behavior 
of the system. Fig. illustrates that shifting the strength 
distributions of beams the functional form of the consti- 
tutive behavior remains the same, however, the value of 
the critical stress and strain vary in a relatively broad 
range. 

The same analysis can also be performed by fixing the 
values Ai and A2 and changing the relative width of the 
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FIG. 8: Fraction of fibers broken by the stretching mode as 
a function of e for different values of the corresponding shape 
parameter mi. Strain controlled simulation with failure due 
to the &R-criterion, N = 90000, Ai = A 2 = 1.0, m 2 = 9. 



two distributions by varying one of the Weibull shape pa- 
rameters m. We find it convenient to shift mi, the shape 
parameter of the stretching mode instead of mi . It can 
be observed in Fig. that for this choice of the scale pa- 
rameters A, the value of the critical strain hardly changes, 
however the critical stress nearly doubles as compared to 
Fig. 

Although the effect on the final fraction of beams bro- 
ken by each mode, see inset of Fig. is not as pronounced 
as for shifting A, we should also consider the fraction of 
fibers broken up to a value of e during the loading process 
(Fig. [SJ. It should be noted that the end points of the 
respective curves in Fig.Elare just the final fraction num- 
bers in Fig. \7\ but the curves show a strong spread for 
intermediate values of e. This demonstrates that chang- 
ing the amount of disorder in the breaking thresholds 
strongly influences the process of damaging of the inter- 
face. 

We apply the methods outlined in the previous para- 
graphs to the von Mises case. Obviously, Eq. (|22|l docs 
not allow for a strict separation of the two modes. How- 
ever, the breaking of a beam at a certain value e c is dom- 
inated by stretching if 

(m* >»m. (33) 

With the previous prescriptions for the failure func- 
tions Eqs. (| 1 3|> - we again find a massive influence on 
the constitutive behavior and the final number of bro- 
ken beams, see Fig. El The inset of Fig. E| demonstrates 
that a crossover between stretching and bending prepon- 
derance occurs also in the von Mises case. 



0.4 




FIG. 9: Constitutive behavior for different values of the bend- 
ing scale parameter A2 . Strain controlled simulation with the 
von Mises criterion , N = 90000, Ai = 1.0, mi = 7712 = 2. 
The inset presents the fraction of beams whose failure was 
dominated by the stretching or bending mode. 



VI. AVALANCHE STATISTICS 

The stress controlled loading of the glued interface is 
accompanied by avalanches of simultaneously failing el- 
ements. The avalanche activity can be characterized by 
the distribution D(A) of burst sizes A defined as the 
number of beam breakings triggered by the failure of a 
single beam. In the framework of simple fiber bundle 
models, it has been shown analytically that global load 
sharing gives rise to a power law distribution of avalanche 
sizes for a very broad class of disorder distributions of 
materials strength [(J 

D(A) oc A' 5 (34) 

with an universal exponent S = 5/2. 

In the previous sections we have shown that in our 
model the interplay of the two breaking modes results in 
a complex failure mechanism on the microscopic level, 
which is strongly affected by the distributions of the 
breaking thresholds. In order to investigate the bursts 
of breaking beams we performed stress controlled simu- 
lations on large systems (N — 10 4 ... 16 • 10 6 ) with both 
the OR and von Mises type breaking criterion. In Figure 
Ell the simulation results are compared to the avalanche 
size distribution of a simple fiber bundle model where 
failure occurs solely due to stretching [E IS El ■ ^ n 
all the cases the avalanche size distributions can be fit- 
ted by a power law over three orders of magnitude. The 
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FIG. 10: Avalanche size distribution D( A) for pure stretching 
of a fiber bundle, and the two beam breaking conditions for 
system sizes N = 16 ■ 10 6 , averaged over 100 runs. A fit 
with the best result D cx A~ 2 ' 56 over almost four decades 
is provided. The inset shows the dependency of the largest 
avalanche A; on the system size for the three cases. Again, 
no difference is found. 



best fit exponent of 5 = 2.56 ± 0.08 was derived from the 
system of size N = 16 • 10 6 beams, with an average taken 
over 100 samples. The size of the largest avalanche in the 
inset of Fig. proved to be proportional to the system 
size. It can be concluded that the beam model belongs 
to the same universality class as the fiber bundle model 

VII. LOCAL LOAD SHARING 

During the failure of interfaces, stress localization is 
known to occur in the vicinity of failed regions, which 
results in correlated growth and coalescence of cracks. 
In our model this effect can be captured by localized in- 
teraction of the interface elements, which naturally oc- 
curs when the two solid blocks are not perfectly rigid 
[l2^ . For simplicity, in our model solely the extremal 
case of very localized interactions is considered, i.e. af- 
ter breaking of a beam in the square lattice, the load is 
redistributed equally on its nearest and next-nearest in- 
tact neighbors. This localized load sharing (LLS) results 
in growing failed regions (cracks) with high stress con- 
centration along their perimeter 0, 0, H3 ■ Figure ^] 
shows the last stable configuration of a beam lattice pre- 
ceding global failure, which was obtained using the OR 
criterion for beam breaking. Due to the stress concentra- 
tion around cracks, the onset of a catastrophic avalanche 
occurs at lower external loads making the macroscopic 
response of the interface more brittle compared to the 
case of global load sharing 0, 0, H3 ■ 

As for global load sharing, we shift the relative im- 
portance of the two breaking modes by changing their 



FIG. 11: Snapshot of a LLS system at the last stable config- 
uration. The color coding represents the load per beam, with 
broken beams carrying a vanishing load. The system size is 
L = 100. 
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FIG. 12: Constitutive curves in the LLS case, shifting A2 and 
keeping the parameters Ai = 1.0 and mi = mi = 2 fixed. 
Stress controlled simulation of N = 10000 fibers averaged 
over 300 runs. 

threshold distributions, and record the influence on 
micro- and macroscopic system properties. We consider 
here the Oi?-critcrion, and use two Weibull distributions 
with parameters Ai, A2 and TOi,m 2 , where the indices 1 
and 2 denote the stretching and bending mode, respec- 
tively. Varying A2 for a fixed Ai, we find a considerable 
influence on the constitutive properties, as Fig. ^| illus- 
trates. 

We investigated also the distribution of avalanche sizes 
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FIG. 13: Distribution of avalanche sizes for LLS for three 
values of A2, with Ai = 1 and mi = m,2 = 2 fixed. Simulations 
were performed using the OR criterion for a bundle of 10000 
beams averaged over 300 runs. 



for LLS, Fig. ^| where we vary the scale parameter A2 
of the bending mode g. We find merely a shifting to 
different amplitudes, but no considerable effect on the 
shape of the distribution function, which is similar to the 
one reported in 0- In comparison to the GLS case, we 
should note that large avalanches cannot occur, and the 
functional form of the curves can be approximated by 
a power law with an exponent higher than for GLS in 
agreement with Refs. |rl Ej . 

VIII. CONCLUDING REMARKS 

Fiber bundle models have been applied to describe var- 
ious aspects of the failure of heterogeneous interfaces. 
However, fibers can sustain solely elongation, and hence 
cannot account for more complex deformation states of 
interface elements, which naturally occurs under shear 
loading. We constructed a novel type of model for the 
shear failure of the glued interface of two solid blocks. In 
the model the interface is discretized in terms of elastic 
beams which experience stretching and bending deforma- 
tion under shear. Breaking of a beam can be caused by 
both deformations resulting in two failure modes. The 
mechanical strength of beam elements is characterized 
by the two threshold values of stretching and bending 
the beam can withstand. The beams are assumed to 
have identical elastic properties, the heterogeneous mi- 
crostructure is represented by the disorder distribution 
of the breaking thresholds. In the model we assume that 
the two solid blocks are perfectly rigid which results in a 



global load redistribution over the intact beams following 
the breaking events. 

We presented a detailed study of the macroscopic re- 
sponse and the progressive damaging of the interface un- 
der quasistatic loading. Making use of the global load 
sharing of intact beams, we obtained the analytic solution 
of the model for the constitutive behavior and the amount 
of damage during the course of loading. In order to ex- 
plore the microscopic process of damaging we worked 
out an efficient simulation technique which enables us to 
study large systems. We demonstrated that the disorder 
distribution and the relative importance of the two failure 
modes have a substantial effect on both the microscopic 
damage process and the macroscopic constitutive behav- 
ior of the interface. Varying its parameters, the model 
provides a broad spectrum of material behaviors. Simu- 
lations showed that the failure of the interface proceeds 
in bursts of simultaneously breaking beams. The distri- 
bution of burst sizes follows power law behavior with an 
exponent equal to the one of simple fiber bundles. Un- 
der stress controlled loading conditions, the macroscopic 
failure of the interface occurs analogously to phase tran- 
sitions, where our beam model proved to be in the same 
universality class as the equal load sharing fiber bundle 
model 0, IM H3 . We showed that the localized inter- 
action of beams leads to a more brittle behavior of the 
interface, which implies a more abrupt transition at the 
critical load. 

Beam models have been successfully applied to study 
the fracture of cohesive frictional materials where cracks 
usually form along the grain-grain interface. Beam el- 
ements proved to give a satisfactory description of the 
interfacial failure of grains and the emerging micro- and 
macro-behavior of materials 25]. Our beam model pre- 
sented here provides a more realistic description of the in- 
terface of macroscopic solid bodies than the simple fiber 
bundle model and is applicable to more complex loading 
situations. Experiments on the shear failure of glued in- 
terfaces are rather limited, especially on the microscopic 
mechanism of the progressive damage, which hinders the 
direct comparison of our theoretical results to experimen- 
tal findings. Our work in this direction is in progress. 
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